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Abstract 



The description of the dynamics at high density parton regime is one of the main 
open questions of the strong interactions theory In this paper we address the shad- 
owing corrections (SC) in the scaling violations of the F2 structure function using the 
eikonal approach. We propose a procedure to estimate the distinct contributions to 
the SC for F 2 and its slope and show that the recent ZEUS data can be described if 
the SC in the quark and gluon sectors are considered. The radius dependence of the 
SC is estimated. Moreover, we calculate the superior limit above which the unitarity 
corrections cannot be disregarded at low Q 2 and show that the recent HERA data 
overcomes this bound. 
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1 Introduction 



The description of the dynamics at high density parton regime is one of the main 
open questions of the strong interactions theory. While in the region of moderate 
Bjorken x (x > 10 -2 ) the well-established methods of operator product expansion and 
renormalization group equations have been applied successfully, the small x region 
still lacks a consistent theoretical framework (For a review see [[[]]). Basically, its is 
questionable the use of the DGLAP equations [Q], which reflects the dynamics at 
moderate x, in the region of small values of x. The traditional procedure of using the 
DGLAP equations to calculate the gluon distribution at small x and large momentum 
transfer Q 2 is by summing the leading powers of a s InQ 2 In(-), where a s is the strong 
coupling constant, known as the double-leading-logarithm approximation (DLLA). 
In axial gauges, these leading double logarithms are generated by ladder diagrams 
in which the emitted gluons have strongly ordered transverse momenta, as well as 
strongly ordered longitudinal momenta. Therefore the DGLAP must breakdown at 
small values of x, firstly because this framework does not account for the contributions 
to the cross section which are leading in a s ln(~) 0. Secondly, because the parton 
densities become large and there is need to develop a high density formulation of 
QCD i. 

There has been intense debate on to which extent non-conventional QCD evolution 
is required by the deep inelastic ep HERA data |], ||. Good fits to the F 2 data for 
Q 2 > 1 GeV 2 can be obtained from distinct approaches, which consider DGLAP 
and/or BFKL evolution equations || [/]]. In particular, the conventional perturbative 
QCD approach is very successful in describing the main features of HERA data and, 
hence, the signal of non-conventional QCD dynamics is hidden or mimicked by a 
strong background of conventional QCD evolution. Our goal in this paper is the 
role of the shadowing corrections (SC) in F 2 and its slope. In the last twenty years, 
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several authors (see || for some phenomenological analysis) have performed a detailed 
study of the shadowing effect although without a strong experimental evidence of 
this effect in the data, mainly since the main observable, the F 2 structure function, 
is inclusive to the effects in the gluon distribution. Recently we have estimated the 
shadowing corrections to the F| and Fl at HERA kinematic region using the eikonal 
approach ||. These observables are directly dependent on the behavior of the gluon 
distribution. We have shown that the shadowing corrections to these observables 
are important, however the experimental errors in these observables are still large to 
allow a discrimination between our predictions and the DGLAP predictions. Here we 
estimate the shadowing corrections to the scaling violations of the proton structure 
function. Basically, there are two possibilities to estimate the SC using the eikonal 
approach. We can calculate damping factors, which represent the ratio between the 
observable with and without shadowing, and subsequently apply these factors in the 
conventional DGLAP predictions. This procedure was used in refs. [TO, IT], also 



considering a two radius model for the nucleon. In this paper we propose a second 
procedure to estimate the SC in DIS, where the observables are directly calculated 
in the eikonal approach and the distinct contributions to the SC are analysed in the 
same approach, reducing the number of free parameters. A larger discussion about 
the distinct procedures is made in section II. 



The recent HERA data on the slope of the F 2 structure function |Iz| present at 
small values of x and Q 2 a different behavior than predicted by the standard DGLAP 
framework. Basically, the HERA data present a 'turn over' of the slope around 
x ~ 10~ 4 , which cannot be described using the GRV94 parametrization jn| and the 



DGLAP evolution equations. We show that this behavior is predicted by the eikonal 
approach considering the shadowing corrections for the gluon and quark sectors. 

The value of the shadowing corrections depends crucially on the size of the target 
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R. The value of the effective radius R depends on how the gluon ladders couple to 
the proton; i.e., on how the gluons are distributed within the proton [ l4| . In this 
paper we estimate the R dependence of the SC. We show that the HERA data on 
the F 2 and its slope can be described consistently using R 2 = 5GeV~ 2 . This value 
agrees with the HERA results on the diffractive J/\l/ photoproduction |15, [16 . 



The steep increase of the gluon distribution predicted by DGLAP and BFKL 



equations at high energies would eventually violate the Froissart bound [17]], which 
restricts the rate of growth of the total cross section to ln 2 (^). This bound may not be 
applicable in the case of particles off-mass shell flTB| , but in this paper we present an 
approach for this problem. Basically, we estimate a limit below which the unitarity 
corrections may be disregarded and show that the recent HERA data surpass this 
boundary, as predicted in |1{J, at small values of x and Q 2 . 

This paper is organized as follows. In section II, the eikonal approach and the 
shadowing corrections for F 2 and its slope are considered. We estimate the distinct 
contributions for the SC and demonstrate that the dF ^ g Q2 - data may be described 
considering the shadowing in the gluon and quark sectors. In section III, we estimate 
the R dependence of the shadowing corrections. In section IV, we present a boundary 
related to unitarity for F 2 and dF ^ g q2 ^ an d show that the actual HERA data for small 
x and Q 2 overcomes this boundary. Therefore, the shadowing corrections should be 
considered in the calculation of the observables in this kinematic region. Finally, in 
section V, we present a summary of our results. 



2 The Shadowing Corrections in pQCD 

The deep inelastic scattering (DIS) is usually described in a frame where the proton 
is going very fast. In this case the shadowing effect is a result of an overlap of the 



5 



parton clouds in the longitudinal direction. Other interpretation of DIS is the intuitive 



view proposed by V. N. Gribov many years ago for the DIS on nuclear targets ||20|| . 
Gribov's assumption is that at small values of x the virtual photon fluctuates into 
a qq pair well before the interaction with the target, and this system interacts with 
the target. This formalism has been established as an useful tool for calculating 
deep inelastic and related diffractive cross section for 7*p scattering in the last years 
pi] , The Gribov factorization follows from the fact that the lifetime of the qq 
fluctuation is much larger than the time of the interactions with partons. According 
to the uncertainty principle, the fluctuation time is ~ where m denotes the target 
mass. 

The space-time picture of the DIS in the target rest frame can be viewed as the 
decay of the virtual photon at high energy (small x) into a quark- ant iquark pair long 
before the interaction with the target. The qq pair subsequently interacts with the 
target. In the small x region, where x -C (R is the size of the target), the qq 
pair crosses the target with fixed transverse distance r t between the quarks. It allows 
to factorize the total cross section between the wave function of the photon and the 
interaction cross section of the quark- ant iquark pair with the target. The photon 
wave function is calculable and the interaction cross section is modeled. Therefore 
we have that the proton structure function is given by 

Q 2 



F 2 (x,Q 2 ) = -f— fdzf d 2 r t \^(z,r t )\ 2 a^(z,r t ) , (1) 



where 



Mz,n)\ 2 = ^f^W + (1 - *)V^iK) 2 +m)Kl{er t f} , (2) 

a em is the electromagnetic coupling constant, e = z(l — z)Q 2 + m 2 , mj is the quark 
mass, nf is the number of active flavors, e 2 is the square of the parton charge (in 
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units of e), K i are the modified Bessel functions and z is the fraction of the photon's 
light-cone momentum carried by one of the quarks of the pair. In the leading log(l/x) 
approximation we can neglect the change of z during the interaction and describe the 
cross section a qq (z,r 2 ) as a function of the variable x. Considering only light quarks 
{i = u, d, s) F 2 can be expressed by [|19j 



u,d,s * 

We have introduced a cutoff in the superior limit of the integration in order to elim- 
inate the long distance (non-perturbative) contribution in our calculations. In this 
paper we assume Ql = OAGeV 2 as in our previous works in this subject. 

We estimated the shadowing corrections considering the eikonal approach [E3], 
which is formulated in the impact parameter space. Here we review the main as- 
sumptions of the eikonal approach. In the impact parameter representation, the 
scattering amplitude A(s, t), where t = —q 2 is the momentum transfer squared, is 
given by 

= 7^ / d 2 q t e-^A(s,t) . (4) 

The total cross section is written as 

a tot (s) = 2 J d\lma(s,b t ) , (5) 

and the unitarity constraint stands as 

2Ima(s,b t ) = \a(s,b t )\ 2 + G in (s,b t ) (6) 

at fixed b t , where Gi n is the sum of all inelastic channels. For high energies the general 
solution of Eq. (10) is: 



a(s,b t ) = i 



"(s,6 t ) 

1 — e 2 



(7) 
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where the opacity Q(s, b t ) is a real arbitrary function, which is modeled in the eikonal 
approach. 

Using the s-channel unitarity constraint (0) in the expression (H), the F 2 structure 



function can be written in the eikonal approach as [24 



1 2 f~k d 2 r t 



Q 2 ) = 7TJ E 4 / r / d\{l - e-^^M } ) (8) 

ZiTT i J rr V+ J 

u,d,s 1 

where the opacity fl qq (x, r t ,b t ) describes the interaction of the qq pair with the target. 

In the region where Q q g is small (Q q q 1) the bt dependence can be factorized 
as flgq = Q q qS(b t ) @] , with the normalization / d 2 b t S(b t ) = 1. The eikonal approach 
assumes that the factorization of the bt dependence Q q g = Q qq S(b t ), which is valid 
in the region where Q qq is small, occurs in the whole kinematical region. The main 
assumption of the eikonal approach in pQCD is the identification of opacity Q qq with 
the gluon distribution. In |19[ the opacity is given by 

Uq-q = ^7c 2 r^xG(x,Q 2 ) , (9) 

where xG(x, Q 2 ) is the gluon distribution. Therefore the behavior of the F 2 structure 

function @ in the small-x region is mainly determined by the behavior of the gluon 

distribution in this region. 

The use of the Gaussian parametrization for the nucleon profile function S(b t ) = 
1 - b2 

where R is a free parameter, simplifies the calculations. In general this 
parameter is identified with the proton radius. However, R is associated with the 
spatial gluon distribution within the proton, which may be smaller than the proton 
radius (see discussion in the next section). 

Using the expression (Q) in (|8|) and doing the integral over b t , the master equation 
for F 2 is obtained p4 



F 2 (x } Q 2 ) = £ e 2 ho -^{C + lniK^r^ + E^K^r 2 ))} , (10) 

u,d,s o 1 1 



where C is the Euler constant, E\ is the exponential function, and the function 
K q (x,r^) = ix r\ xG(x, p). Expanding the equation (|10|) for small K g , the first 
term (Born term) will correspond to the usual DGLAP equation in the small x region, 
while the other terms will take into account the shadowing corrections. 

The slope of F 2 structure function in the eikonal approach is straightforward from 



the expression (|10|) . We obtain that 



dF ^ g Qp = E ej{C + ln( Kq (xy t )) + E 1 ( Kq (xy t ))} . (11) 



u,d,s 



The expressions (0) and (|H|) predict the behavior of the shadowing corrections 
to F2 and its slope considering the eikonal approach for the interaction of the qq with 
the target. In this case we are calculating the SC associated with the passage of the 



qq pair through the target. Following [ljj we will denote this contribution as the 
quark sector contribution to the SC. 

The behavior of F 2 and its slope are associated with the behavior of the gluon 
distribution used as input in (10) and (|TT1). In general, it is assumed that the gluon 



distribution is described by a parametrization of the parton distributions (for example: 
GRV, MRS, CTEQ) |25|, p6fl . In this case the shadowing in the gluon distribution 
is not included explicitly. In a general case we must also estimate the shadowing 
corrections for the gluon distribution, i.e. in the quark and the gluon sectors. In this 
case we must estimate the SC for the gluon distribution using the eikonal approach, 
similarly to the F 2 case. This was made in [57[ and here we only present the main 
steps of the approach. 

The gluon distribution can be obtained in the target rest frame considering the 
decay of a virtual gluon at high energy (small x) into a gluon-gluon pair long before 
the interaction with the target. The gg pair subsequently interacts with the target, 
with the transverse distance r t between the gluons assumed fixed. In this case the 
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cross section of the absorption of a gluon g* with virtuality Q 2 can be written as 

^nudcon^) = f l dz f fn^g* (Q2 ;ri)X)Z) |2 a93+ nucleon ( ^ r2) (12) 

Jo J n 

where z is the fraction of energy carried by the gluon and is the wave function of 
the transverse polarized gluon in the virtual probe. Furthermore, a 99+nuclcon (z, r 2 ) is 
the cross section of the interaction of the gg pair with the nucleon. Considering the 
s-channel unitarity and the eikonal model, equation (Jl^) can be written as 

a g*+™i™^Q*) = J 1 dz jfll J^.\*f((f,r t ,x,z)\ 2 (l-e-I^W) , 

where the factorization of the b t dependence in the opacity Q gg (x, r t , b t ) was assumed. 
Using the relation a 9 * +nucleon (x, Q 2 ) = 47 ^" a xG(x, Q 2 ) and the expression of the wave 
^/ 9 calculated in PSL |2"7fl , the Glauber-Mueller formula for the gluon distribution is 
obtained as 



2 4 f^dx' [°°d 2 r t f°°d 2 b t 
xG{x,Q ) = — — — 2 

7V Z Jx X' J-^r TTVf Jo 7T 



1 - e -2°N 



la^(x'4)S(bt) 



(13) 



where Q gg = ojy describes the interaction of the gg pair with the target. Using the 
Gaussian parametrization for the nucleon profile function, doing the integral over bt, 
the master equation for the gluon distribution is obtained as 

9T? 2 fl tiv* r-^r d 2 r+ 

xG(x,Q 2 ) = ^ ^ / Q o a -^{C + ln(K G (x',r 2 )) + E 1 (k g (x' , r 2 ))} , (14) 
it z Jx x' J-L Tirf 

where the function k g (x, r 2 ) = 7rrf xG(x, ^). Again, if equation ( |i~4"| ) is expanded 
for small k g , the first term (Born term) will correspond to the usual DGLAP equation 
in the small x region, while the other terms will take into account the shadowing 
corrections. The expressions (|10|) , ([TT| ) and (14]) are correct in the double leading 
logarithmic approximation (DLL A). As shown in [53] the DLLA does not work quite 
well in the accessible kinematic region (Q 2 > OAGeV 2 and x > 10~ 6 ). Consequently, 
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a more realistic approach must be considered to calculate the observables. In |23 
the subtraction of the Born term and the addition of the GRV parametrization were 
proposed to the F 2 and xG cases. In these cases we have 



F 2 {x,Q 2 ) = F 2 {x,Q 2 )[Eq. (10)] -F 2 (x,Q 2 )[Born] +F 2 (x,Q 2 )[GRV] , (15) 

and 

xG(x, Q 2 ) = xG(x, Q 2 ) [Eq. (14)] - xG(x, Q 2 ) [Born] + xG(x, Q 2 ) [GRV] , (16) 

where the Born term is the first term in the expansion in n q and n g of the equations 
( |TD| ) and (|T1"D, respectively (see || for more details). Here we present this procedure 



for the F 2 slope. In this case 
dF 2 (x,Q 2 ) dF 2 (x,Q 2 ) ^ dF 2 (x,Q 2 ) ^ 1 dF 2 (x,Q 2 ) 

= ^g^ [Eq - - ^%o^ [Born] + ^q^ [grv] ' (17) 

where the Born term is the first term in the expansion in n q of the equation (|il"|). 
The last term is associated with the traditional DGLAP framework, which at small 
values of x predicts 

dF 2 (x,Q 2 ) 10a s (Q 2 ) [i-* x ( x 2 \ 

-H^ = ^^L dzP q9 (z)—g(—,Q), (18) 

where a s (Q 2 ) is the running coupling constant and the splitting function P qg (x) gives 
the probability to find a quark with momentum fraction x inside a gluon. This 
equation describes the scaling violations of the proton structure function in terms of 
the gluon distribution. We use the GRV parametrization as input in the expression 
(|18). In the general approach proposed in this paper we will use the solution of the 
equation ([TBI) as input in the first terms of (]15D and (|17D . As the expression (|i~6|) 
estimates the gluon shadowing, the use of this distribution in the expressions (O) 



and (|TTD, which consider the contribution to SC associated with the passage of qq 
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pair through the target, allows to estimate the SC to both sectors (quark + gluon) 
of the observables. Our goal is the discrimination of the distinct contributions to the 



SC in F 2 and d -^gl. 



In Fig. D we present our results for the F 2 structure function as a function of 
the variable In(-) for different virtualities. We have used R 2 = 5GeV~ 2 in these 
calculations. In the next section the R dependence of our results is analysed. We 
present our results using the expression ( p3|) (quark sector) and using the solution 
of the equation fll6|) as input in the first term of (0) (quark + gluon sector). The 
predictions of the GRV parametrization are also shown. We consider the HERA 
data at low Q 2 since for Q 2 > 6 GeV 2 the SC start to fall down (For a discussion 
of the SC to F 2 considering the quark sector see p|, p9|). We can see that at small 
values of Q 2 the predictions for F 2 considering the quark and the quark-gluon sector 
are approximately identical. However, for larger values of Q 2 the predictions of the 
quark-gluon sector disagree with the HI data fl30| . Therefore, the contribution of 
the gluon shadowing to F 2 in an eikonal approach superestimates the shadowing 
corrections at large Q 2 values. 

In Fig. ^| we present our results for the SC in the dF ^Q2 ^ as a function of x. The 
ZEUS data points [12] correspond to different x and Q 2 value. The (x, Q 2 ) points are 
averaged values obtained from each of the experimental data distribution bins. Only 
the data points with < Q 2 > > 0.52 GeV 2 and x < 10 _1 were used here. 

The SC are estimated considering the expression QT7D (quark sector) and using 
the solution of the equation (|16D as input in the first term of (|T7|) (quark + gluon 
sector). Moreover, the predictions of the traditional DGLAP framework, which at 
small values of x is given by the expression (|I8|) are also presented. We can see 
that the DGLAP predictions fail to describe the ZEUS data at small values of x and 
Q 2 . However we see that in the traditional framework (DGLAP + GRV94) a 'turn 
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over' is also present at small values of x and Q 2 . Basically, this occurs since the 
smaller Q 2 value used (< Q 2 >= 0.52 GeV 2 ) is very near the initial virtuality of the 
GRV parametrization, where the gluon distribution is 'valence like'. Therefore the 
gluon distribution and the F 2 slope are approximately flat in this region. For the 
second smaller value of Q 2 used (< Q 2 >= 1.1 GeV 2 ) the evolution length is larger, 
which implies that the gluon distribution (and the F 2 slope) already presents a steep 
behavior. The link between these points implies the 'turn over' presented in Fig. 0. 
The main problem is that this 'turn over' is higher than observed in the ZEUS data. 
This implies that xG(x, Q 2 ) differs from the previous standard expectations in the 
limit of small x and Q 2 . This effect is not observed in the F 2 structure function since 
it is inclusive to the behavior of the gluon distribution, which can be verified analysing 
the predictions of the distinct parametrizations. The gluon distribution predicted by 
these parametrizations differs in approximately 50 %. 

The prediction of the gluon sector, which is obtained using the solution of the 
expression (|T6| ) as input in ([TJ) is also presented. We can see that at larger values of 
Q 2 and x all predictions are approximately identical. However, at small values of x 
and Q 2 , the ZEUS data is not well described considering only the quark or the gluon 
sector to the SC. The contribution of the gluon shadowing is essential in the region of 
small values of x and Q 2 , i.e. a shadowed gluon distribution should be used as input 
in the eikonalized expression (|T7|) in this kinematic region. 

Our conclusion is that at small values of x and Q 2 it should be considered the 
contribution of the gluon shadowing to estimate the SC to F 2 and its slope in the 
eikonal approach. While for F 2 the contribution of the gluon shadowing may be 
disregarded, it is essential for the F 2 slope. The dF ^ g Qi ^ data show that a consistent 
approach should consider both contributions at small x and Q 2 . 

Before we conclude this section some comments are in order. We show that the 
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dF (UogQ 2 " data can be successfully described considering the shadowing corrections in 



the quark and gluon sectors. A similar conclusion was obtained in [10, 11], where the 
eikonal approach was also used to estimate the SC in the quark and gluon sectors, 
but a distinct procedure was used to estimate the SC for the F 2 slope. In [[TI|] damp- 
ing factors are calculated separately for both sectors and applied to the standard 
DGLAP predictions. The behavior of the gluon distribution at small values of Q 2 
was modeled separately, since the gluon distribution ( 14 ) vanish for Q 2 = Q\- This 
procedure introduces a free parameter /i 2 , beyond the usual ones used in the eikonal 
approach (Qq, R 2 ). The distinct procedure proposed here estimates the observables 
directly within the eikonal approach and the shadowing corrections in the different 
sectors are calculated within the same approach. In our calculations there are only 
two free parameters: (i) the cutoff (Qq = 0.4 GeV 2 ) in order to eliminate the long 
distance contribution, and (ii) the radius R (R 2 = 5GeV~ 2 ). The choice of these 
parameters is associated with the initial virtuality of the GRV parametrization used 
in our calculations, and the estimates obtained using the HERA data on diffractive 
photoproduction of J/ty vector meson (see discussion in the next section) respec- 
tively |TR |I6|| . In our procedure the region of small values of Q 2 ~ Qq is determined 
by the behavior of the GRV parameterization in this region, since we are using the 
eq. ( |1"6D to calculate the gluon distribution. For Q 2 = Ql the two first terms of 
(|16|) vanish and the gluon distribution is described by the GRV parameterization, i. e. 
xG{x,Q 2 ) = xG{x,Q 2 )[GRV}. 

The eikonal approach describes the ZEUS data, as well as the DGLAP evo- 
lution equations using modified parton distributions. Recently, the MRST group 



[ 31 has proposed a different set of parton parametrizations which consider a ini- 
tial 'valence-like' gluon distribution. This parametrization allows to describe the F 2 
slope data without an unconventional effect. This occurs because there is a large 
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freedom in the initial parton distributions and the initial virtuality used in these 
parametrizations. We believe that only a comprehensive analysis of distinct ob- 
servables (Fl, F$, dF ^ g Q2 ~ ) will allow a more careful evaluation of the shadowing 
corrections at small x P, p9| . 



3 The radius dependence of the shadowing cor- 
rections 



The value of SC crucially depends on the size of the target |jT4|| . In pQCD the value 
of R is associated with the coupling of the gluon ladders with the proton, or to put 
it in another way, on how the gluons are distributed within the proton. R may be 
of the order of the proton radius if the gluons are distributed uniformly in the whole 
proton disc or much smaller if the gluons are concentrated, i.e. if the gluons in the 
proton are confined in a disc with smaller radius than the size of the proton. 

Considering the expression (||), assuming Q q g < 1 and expanding the expression 
to 0(Vl 2 ) we obtain 

^.« 2 )^E4/Jf /^{in*-^} . (i9) 

u,d,s q2 t 

Using the factorization of the opacity and the normalization of the profile function 
we can write F 2 as 

FiM 2 ) = A E e / /5 7? {¥^-\^S d 2 btS 2 (b t )\ . (20) 

u,d,s ~qZ * 

The second term of the above equation represents the first shadowing corrections 
for the F 2 structure function. Assuming a Gaussian parametrization for the profile 
function we obtain that the screening is inversely proportional to the radius. Therefore 
the shadowing corrections are strongly associated with the distributions of the gluons 
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within the proton. In this section we estimate the radius dependence of the shadowing 
corrections, considering the F 2 and dF ^ g q2 ~ data. First we explain why the radius is 
expected to be smaller than the proton radius. 

Consider the first order contribution to the shadowing corrections, where two lad- 
ders couple to the proton. The ladders may be attached to different constituents of 
the proton or to the same constituent. In the first case the shadowing corrections are 
controlled by the proton radius, while in the second case these corrections are con- 
trolled by the constituent radius, which is smaller than the proton radius. Therefore, 
on the average, we expect that the radius will be smaller than the proton radius. 
Theoretically, R 2 reflects the integration over bt in the first diagrams for the SC. 

In Fig. |H we present the ratio 

F 2 (:r,Q 2 )[Eq. (15)] 
F 2 (x,Q 2 )[GRV] ' 1 } 

where F 2 (x, Q 2 )[GRV] = E u ,d,s e } N(^<? 2 ) + xq(x,Q 2 )] + F^(x,Q 2 ) is calculated 
using the GRV parametrization. For the treatment of the charm component of the 



structure function we consider the charm production via boson-gluon fusion [13]. In 
this paper we assume m c = 1.5 GeV. In Fig. |] we present the ratio 



^IgPlGRV] • (22) 

The function ^ffl [GRV] was calculated using the expression fll8f) and the GRV 



parametrization. Our results are presented as a function of ln(-) at different virtu- 
alities. We can see that the SC are larger in the ratio R s and that our predictions of 
SC are strongly dependent of the radius R. Moreover, we see clearly the SC behavior 
inversely proportional with the radius. 

In Fig. [| we compare our predictions for the SC in the F 2 structure function and 
the HI data |3(| as a function of In(-) at different virtualities and some values of the 
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radius. Our goal is not a best fit of the radius, but eliminate some values of radius 
comparing the predictions of the eikonal approach and HERA data. We consider 
only the quark sector in the calculation of SC, which is a good approximation in 
this observable, as shown in the previous section. The choice R 2 = 1.5 GeV 2 does 
not describe the data, i.e. the data discard the possibility of very large SC in the 
HERA kinematic region. However, there are still two possibilities for the radius which 
reasonably describe the F 2 data. To discriminate between these possibilities we must 
consider the behavior of the F 2 slope. 

In Fig. H we present our results for dF ^ x g Q 2 ^ considering the SC only in the quark 
sector. Although in the previous section we have demonstrate that the contributions 
of the quark and gluon sectors should be considered, here we will test other possibili- 
ties to describe the data: the dependence on the radius R. Our results show that the 
best fit of the data occurs at small values of R 2 , which are discarded by the F 2 data. 
Therefore, in agreement with our previous conclusions, we must consider a general 
approach to describe consistently the F 2 and dF ^ g q2 ^ data. In Fig. [7] we present 
our results for dF ^ g Qi - considering the SC in the gluon and quark sector for different 
values of R 2 , calculated using the general approach proposed in the previous section. 
The best result occurs for R 2 = 5GeV~ 2 , which also describes the F 2 data. 

The value for the squared radius R 2 = 5 GeV~ 2 obtained in our analysis agrees 
with the estimates obtained using the HERA data on diffractive photoproduction of 
J/\l/ meson [F5, IE]. Indeed, the experimental values for the slope are B e i = 4GeV~ 2 
and Bi n = 1-66 GeV~ 2 and the cross section for J/^> diffractive production with and 
without photon dissociation are equal. Neglecting the t dependence of the pomeron- 
vector meson coupling the value of R 2 can be estimated [|19|]. It turns out that 
R 2 ~ 5GeV~ 2 , i.e., approximately 2 times smaller than the radius of the proton. 

As an additional comment let us say that the SC to F 2 and its slope may also be 
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analysed using a two radii model for the proton |TTJ. This analysis is motivated by 
the large difference between the measured slopes in elastic and inelastic diffractive 
leptoproduction of vector mesons in DIS. An analysis using the two radii model for 
the proton is not a goal of this paper, since a definite conclusion on the correct model 
is still under debate. 

The summary of this point is that the analysis of the F 2 and dF ^ g Q2 ^ data using 
the eikonal model implies that the gluons are not distributed uniformly in the whole 
proton disc, but behave as concentrated in smaller regions. This conclusion motivates 
an analysis of the jet production, which probes smaller regions within the proton, 
using an approach which considers the shadowing corrections. 



4 A screnning boundary 

The common feature of the BFKL and DGLAP equations is the steep increase of the 
cross sections as x decreases. This steep increase cannot persist down to arbitrary 
low values of x since it violates a fundamental principle of quantum theory, i.e. the 
unitarity. In the context of relativistic quantum field theory of the strong interactions, 
unitarity implies the cross section of a hadronic scattering reaction cannot increase 
with increasing energy s above log 2 s: the Froissart's theorem [ I7| . The Froissart 



bound cannot be proven for off- mass-shell amplitudes [[HJ , which is the case for deep 
inelastic scattering f32fl . Our goal in this section is by using the s-channel unitarity 
(P) and the eikonal approach to estimate a superior limit from which the shadowing 
corrections cannot be disregarded in F 2 and its slope. 

Considering the expression @ for the F 2 structure function, we can write a b t 
dependent structure function given by 

F 2 (x,Q 2 ,b t ) = ^E4 e-in-^W} . (23 ) 



u,d,s 



The relation between the opacity and the gluon distribution (^) obtained in [|T9]], 
is valid in the kinematical region where Q <^ 1. In the eikonal approach for pQCD we 
make the assumption that the relation (§) is valid in all kinematic region. To obtain 
an estimate of the region where the SC are important we consider a superior limit 
for the expression (^), which occurs for Q ^> 1. In this limit the second term in the 
above equation can be disregarded. As the shadowing terms are negative and reduce 
the growth of the F 2 structure function, disregarding the shadowing terms we are 
estimating a superior limit for the region where these terms are not important, i.e. 
a screnning boundary which establishes the region where the shadowing corrections 
are required to calculate the observables. 

The b t dependent structure function in the limit Q 2> 1 is such that 

F^,Q\b t )<^Y.^jt° d 4- (24) 

u,d,s t 

Making the assumption that the bt dependence of the structure function is factorized 
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F 2 (x,Q 2 ,b t )=F 2 (x,Q 2 )S(b 



t) , 



and considering a Gaussian parametrization for the profile function and its value for 
b t = we get (n/ = 3) 



As a result 



p2 

2\ ^ - rt (r\1 /02\ 



F 2 (x,Q 2 ) < —(Q-Qo 



3tt 2 
R 2 Q 2 

3tt 2 



< (26) 
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The above limit is our estimate for the screnning boundary for the F 2 structure 
function. 

The screnning boundary for the F 2 slope is straightforward from the expression 
(H). We get 

dF 2 (x,Q 2 ) W 
dlogQ 2 3tt 2 ' 1 ' 



This expression agrees with the expression obtained in |19 



Clearly expressions (|2q ) and (|27| ) serve only as a rough prescription for estimating 
the region where the corrections required by unitarity cannot be disregarded. A more 
rigorous treatment would be desirable, but remains to be developed. 

Using the above expressions we can make an analysis of HERA data. We use 
R 2 = 5 GeV~ 2 in the calculations. In Fig. || we compare our predictions with the 
F2 data from the HI collaboration. We can see that data at larger values of Q 2 



(Q 2 > 8.5 GeV 2 ) do not violate the limit (|26|) . However, the data at smaller values 
of Q 2 and x violate this limit. This indicates that we should consider the SC for this 
kinematical region. In Fig. [| we present our results for the F2 slope. We see that the 
data for small Q 2 and x (Q 2 < 2.5 GeV 2 , x < 10~ 4 ) violate the limit fl27|) , stressing 
the need of the shadowing corrections. Therefore for small values of x and Q 2 the 
observables must be calculated using an approach which takes them into account. 



5 Summary 

In this paper we have presented our analysis of the shadowing corrections in the 
scaling violations using the eikonal approach. We shown that the dF ^ g Q2 ^ data can 
be described successfully considering the shadowing corrections in the quark and gluon 
sectors. Furthermore, we have considered the radius dependence of these corrections 
and an unitarity boundary. From the analysis of the R dependence of the SC in the 



20 



eikonal approach we have shown that the value R 2 = 5 GeV~ 2 allows to describe 
the HERA data. This value agrees with the estimate obtained independently in the 
diffractive J/^f photoproduction. Using the eikonal approach and the assumption of 
b t factorization of the F 2 structure function a screnning boundary is analysed. This 
boundary constrains the region where the corrections required by unitarity may be 
disregarded; or in other words, a limit for applicablity of standard perturbative QCD 
framework. We have shown that the HERA data at small x and Q 2 violate this limit, 
which implies that the shadowing corrections are important in the HERA kinematic 
region. 

Our conclusion is that the shadowing effect is important already at HERA kine- 
matic region. We believe that the analysis of distinct observables (Fx,, F 2 C , dF ^ x g Q 2 ^ ) 
at small values of x and Q 2 will allow to evidentiate the shadowing corrections. 
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Figure 1: The structure function as a function of the variable ln(-) for different 



virtualities. Data from HI [30]. The solid curve corresponds to GRV, the dashed 
curve to SC in the quark sector, the long-short dashed curve to SC in both quark and 
gluon sectors. 
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Figure 2: The F2 slope as a function of the variable x. Data from ZEUS [[Cg] . The 
data points correspond to a different x and Q 2 value. The solid curve corresponds to 
SC in both quark and gluon sectors. 
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Figure 3: The ratio Ro = F2<yX '® as a function of the variable In(-) for 

different virtualities and radii. 
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Figure 4: The ratio i? s = Q H) rr , DUl — as a function of the variable In(^) for 



different virtualities and radii. 
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Figure 5: The F 2 structure function as a function of the variable In(-) for differ- 
ent virtualities and radii. Only the shadowing corrections in the quark sector are 
considered. Data from HI |3(J . 
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Figure 6: The F2 slope as a function of the variable x for different radii. Only the 
shadowing corrections in the quark sector are considered. Data from ZEUS ||12|| . The 
data points correspond to a different x and Q 2 value. 
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X 



Figure 7: The F2 slope as a function of the variable x for different radii. The shadow- 



ing corrections in the gluon-quark sector are considered. Data from ZEUS fll2| . The 
data points correspond to a different x and Q 2 value. 
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Figure 8: Screnning boundary to the F2 structure function. For points above the 
boundary the corrections cannot be disregarded. Data from HI . 
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